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ABSTRACT 

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is 
an important tracer of cluster merger activity, and is the result of complex physical processes that in- 
volve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known 
as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a 
magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this 
emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution 
of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of 
cosmic ray electron acceleration that are then input into radio emission models. We have determined 
the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio obser- 
vations to compare to observed galaxy clusters. We find a significant dependence of the observed 
morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding 
the interpretation of observed radio features in clusters. We also find that a given shock should not 
be characterized by a single Mach number. We find that the bulk of the radio emission comes from 
gas with T > 5 x 10 7 K, p ~ 10~ 28 — 10~ 27 g/cm 3 , with magnetic field strengths of 0.1 — 1.0/zG and 
shock Mach numbers of A4 ~ 3 — 6. We present an analysis of the radio spectral index which suggests 
that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the 
polarization fraction and position angle of the simulated radio features, and compare to observations. 
Subject headings: cosmology: theory — magnetohydrodynamics — methods: numerical — cosmic 
rays — radiation mechanisms: nonthermal 



1. INTRODUCTION 

Galaxy clusters are hosts to a variety of thermal and 
non-thermal phenomena, many of which are the result 
of cosmological structure formation. The study of rela- 
tivistic particles in galaxy cluster environments was mo- 
tivated by the observat i on of the radio halo in the Coma 
cluster by Large et al. (1959), and has since grown into 
an industry ot observations, theory, and simulation. For 
a r eview on current radio observations o f galax y clusters 
see |Ferrari et all (|2008|; iFeretti et all fl2012|), and for 



a revi ew on the non-thermal processes see Dolag et al 



(120081). Here we review the basic characteristics of galaxy 
cluster radio "h alos" and giant radio "relics," to use the 
classification in Ferrari et al.| ([2008). Radio halos are 
usually ~ Mpc— sized features in galaxy clusters, closely 
following the X-ray morphology in the central regions 
of the cluster. They are generally characterized by very 
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low (< few percent) linear polarization fractions, and arc 
found in galaxy clusters wit h disturbed morphology and 
no evidence for a, cool core (|Ferrari et al .||2008[ IGiovan- 



nini et al.||2009[ |Feretti et al.||2012| . The origin of the 
emission is thought to be from relativistic (7 ~ 10 4 ) elec- 
trons emitting synchrotron radiation. The source of the 
energy in these electrons, however, is debated. It may 
originate from the decay of pions (the "secondary" or 
"hadronic" model), created by interactions be tween cos- 
mic r ay protons and t he th ermal population (jDennison 
1980| |Dolag fc Enfilin|[2^00] |Miniati et aLl|200l[ ), which 



would be strengthened by the observation ot gamma- 
ray emission in cluster cores. However, initial studies 



of many galaxy clust ers using the FERMI satellite ( Ack 



ermann et al. 2010), as well as for fewer objects with 
other instruments (e.g. MAGIC observations of Per seus 
lAleksic et al.|2010[|2012[ ) , combined with radio data pel 
Itema & Frofumo 2011 



Brunetti et al. 
the energy in cosmic rays to be very 



2012) constrain 
mTjk 10%) of 



the thermal energy in most cases. Others believe that 
the electrons arc turbulcntly accelerated either from the 
thermal population or from aging populations of elec- 
trons eith er from shock acceleration or AGN/supernova 
injection (Brunetti & Lazarian 20111. 

Radio relics, on the other hand, are thought to be ac- 
celerated by first-order Fermi acceleration thro ugh the 
process of Diffusive Shock Acceleration (DSA) ( |Bland 



ford & Ostriker 1978 ). They have a relatively steep radio, 
and therefore inferred electron, spectrum where S cx u~ a 
with a w 1 — 2. These radio sources are not associated 
with any of the cluster galaxies or AGN bubbles. They 
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are also not associated with any point sources in other 
wavelengths, and are usually found in the outskirts of 
clusters. Their location can be up to ~ 2 Mpc from the 
cluster core, and can be extended up t o ~ 1 .5 Mpc in 
length (Ivan Weeren et al.||2010l |2011cl 12012b. In some 



cases these radio features are coincident with X-ray sur 
face brightness and temperature jump s, potentially indi 
eating the presence of a sho ck front (Finoguenov et al. 



2010 



Akamatsu et al. 2012 1. While more rare, double 



radio relics are observed in several systems. Double ra- 
dio relics are unique in that they provide tighter con- 
straints on the geometry and kin ematics of the merging 
clusters( "van Weeren et al.||2011a ). Upcoming radio tele- 
scopes such as LOFAR, the Jansky VLA, and eventually 
the SKA will provide an increase in sensitivity and res- 
olution (both spectral and spatial) that will allow for 
discoveries in blind surveys. Because of this, we are at 
an important time to use simulation and theory to pre- 
dict the number and the properties of relics in cosmo- 
logical samples. Past simulations have focused on both 



single clusters (Rocttiger et al. 1999 Pfrommer et al. 



ters (Hocft et al. 2008 



JSkillman et al.|2011 


Vazza et al. 


21). Both are needed 


in order to 



constrain the plasma physics and how varying environ- 
ments lead to observational quantities such as luminosity 
functions. 

In this paper we investigate the origins, properties, and 
observational implications of a merging galaxy cluster us- 
ing a numerical simulation. For the first time, we start 
from cosmological initial conditions and self-consistently 
evolve the cluster magnetic field from an AGN source the 
equations of magnetohydrodynamics rather than assum- 
ing a magnetic field strength and topology. This allows 
us to explore one scenario in which the magnetic field 
forms, evolves, and interacts with the radio relic emis- 
sion. We describe, in detail, the plasma environment of 
the radio-emitting regions. 

After investigating the properties of the cluster gas, 
we analyzed the resulting radio emission using novel ap- 
proaches to explore systematic effects present in current 
radio observations. We used a new tool to view this Eu- 
lerian grid simulation from arbitrary directions in order 
to demonstrate the effect of viewing angle on the derived 
properties. We then developed the capability to inte- 
grate the polarized radio emission along the line of sight 
to provide the closest comparison to observations. We 
then use this technique to produce polarization fraction 
and position angle maps from our MHD AMR simula- 
tion, and provide comments on the relevance of our re- 
sults to current observations of radio features in galaxy 
clusters. Finally, we discuss the impact of using previous 
assumptions about the magnetic field compared to the 
values that are self-consistently evolved from an AGN 
source. We use this to provide insight into observational 
results. 

2. METHODS 

2.1. Simulations 

Our simulation was r un using a modified version of the 
Enzo cosmology co de ( Bryan k Nor man 199 7a|b Nor- 
man k Bryan|1999||0^hea et al.|2004| ). Enzovi sesblocE 
structured adaptive mesh refinement (AMR; |Berger k 



Colella 1989 1 as a base upon which it couples an Eulerian 
hydrodynaniic solver f or the gas with an N-Body parti 



cle mesh (PM) so lver (Efstathiou et al. 
k Eastwood||1988 ) for the dark ma tter 



1985 Hockney 



n this work we 



utilize the MHD solver described in Collins et al. (2010| 
The solver emplo yed here is spatially second o rder, whi. 



the PPM solver ( |Colella k Woodward! 1984[ ) commonly 
used in Enzo is spatially 3rd order. The net effect on 
the shock-finding algorithm will be to broaden a single 
shock by a small amount. However it will be impossi- 
ble to disentangle this effect from the changes in shock 
structure due to the addition of magnetic forces in the 
evolution. None of the results we present here will be 
sensitive to these small differences. We have extended 
this version of Enzo to inclu de temperature-jump b ased 
shock-fin ding as described in Skillman et al. ( 2008 1 and 
used in |Skillman et al. (2011 ). 

The galaxy cl uster studied in t his work is the same 
as cluster Ul in |Xu et al. (2011). In this work, clus- 



ters were formed from cosmological initial conditions 
and magnetic fields were injected by the most massive 
galaxy at a variety of stages in the cluster evolution. It 
was found that different injection parameters of magnetic 
fields have little impact on the cluster formation history. 
This simulation models the evolution of dark matter, 
baryonic matter, and magnetic fields self-consistently. 
The simulation uses an adiabatic equation of state for 
gas, with the ratio of specific heat being 5/3, and does 
not include heating or cooling physics or chemical re- 
actions. While studies have been done including these 
physical models and their role in charact erizing shocks 



(Kang et al. 2007 Pfrommer et al. 2007), we chose to 



ignore them due to both computational cost as well as 
possible confusion between structure formation shocks 
and those arising from star/galaxy feedback. Addition- 
ally, Kang et al. (2007) found little effect on the overall 
kinetic energy dissipation between simulations with adi- 
abatic gas physi cs and those including cooling and feed- 
back, and while |Pfrommer et al. ( 2007 1 show changes at 



high Mach number, as we will see these have little con 
sequence for the shocks involved with producing radio 
relics. 

The initial conditions o f the simulation are gen erated 



(1999) power 



at redshift z = 30 from an Eisenstein k Hu 
spectrum of density fluctuations in a ACDM universe 
with parameters h = 0.73, Q m = 0.27, fi 6 = 0.044, 
Oa = 0.73, as — 0.77, and n s = 0.96. These parame- 
ters are close to th e values from WMAP3 observations 
( Spergel et al.|2003 |. While these parameters differ from 
the latest constraints, it is largely irrelevant for this par- 
ticular project. The simulated volume is (256 /i _1 Mpc) 3 , 
and it uses a 128 3 root grid and 2 nested static grids 
in the Lagrangian region where the cluster forms. This 
gives an effective root grid resolution of 512 3 cells (~ 
0.69 Mpc) and dark matter particle mass resolution of 
1.07 x lO lo Af . During the course of the simulation, 
8 levels of refinements are allowed beyond the root grid, 
for a maximum spatial resolution of 7.8125 hr x kpc. The 
AMR is applied only in a region of (~ 43 Mpc) 3 where 
the galaxy cluster forms near the center of the simula- 
tion domain. The AMR criter ia in this simulation are 
the same as in Xu et al. (20111. During the cluster for- 
mation but before the magnetic fields are injected, the 
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refinement is only controlled by baryon and dark matter 
density, refining on overdensities of 8 for each additional 
level. After magnetic field injections, in addition to the 
density refinement, all the regions where magnetic field 
strengths are higher than 5 x 10 _8 G are refined to the 
highest level. The importance of using this magnetic 
field refinement criterion in cluster MHD simulations is 



discussed in Xu et al 



(2010) 



The magn etic held initializatio n used is the same 



method in Xu et al. (2008 



2009) as th e origi nal mag- 



netic tower model proposed By" Ei et al. 



J2006L and as- 
sumes the magnetic fields are from the outburst of AGN. 
The magnetic fields are injected at redshift z — 3 in two 
proto-clusters, which belong to two sub-clusters. The in- 
jection locations are the same lo cations in simulations 



Ula and Ulb in |Xu et al | ( |2011| ). There is - 6 x 10 59 
erg of magnetic energy placed into the ICM from each 
injection, assuming that ~ 1 percent of the AGN out- 
burst energy of a sever al 10 s Mr SM BH is in magnetic 
fields. Previous studies ( Xu et al.|2010 ) have shown that 
the injection redshifts and magnetic energy have limited 
impact on the distributions of the ICM magnetic fields 
at low redshifts. 

The simulated cluster is a massive cluster with its basic 
properties at redshift z = as follows: Rviriai = 2.5 
Mpc, 
2.7 x 



% i (total) = 1.9 x 10 15 Mq, M OT „ a/ (gas) 



M„ 

10 14 M w , and T 



virial 



10.3 keV. 



This cluster 

is in an unrelaxed dynamical state at z = with its 
two magnetized sub-clusters of similar size undergoing a 
merger. The total magnetic energy in the simulation at 
z = is 9.6 x 10 60 erg, nearly all of which is within 
the cluster virial radius. The details about the cluster 



formation are described in Xu et al. (2011 1. 



2.2. 



Synchrotron Emission 
We use the sam e technique as was presented i n[Skill- 



man et al. (2011 ) and based on Hoeft & Briiggen 
except we no longer rely on the assumption that t 



pooTf ; 

ne mag- 



netic field is a simple function of density and instead use 
the magnetic field from the simulation. This method as- 
sumes that a fraction of the incoming kinetic energy of 
the gas is accelerated by the shock up to a power-law dis- 
tribution in energy, which extends from the thermal dis- 
tribution. This distribution is that predicted by di ffusive 
shock acceleration theory in the test-particle l imit (f Drury 
19831 |Bell||1978l IMalkov fc O'C Drury|[200T] lAchterberg 
fc Wiersma|2007l | Blandford &: Ostriker|l978lfo'C Drury] 



2012 1 . At the high energy end, it also assumes that there 



is an exponential cutoff determined by the balance of ac- 
celeration and cooling. The total radio power from a 
shock wave of area A, frequency v b s , magnetic field B, 
electron acceleration efficiency £ e , electron power-law in- 
dex s (n e oc E~ s ) l post-shock electron d ensity n e and 
temperature T 2 is ( |Hoeft & Bruggen||2007 ) 

dP{v b 



dv 



6.4 x 10 34 er 5 s" 1 Hz' 1 -. 



A 



Mpc 2 10- 4 cm- 3 

£e / v obs \-s/2 

OMAAGHz' 

{B/iiG) 1+ WV 



x ( T2 ) 3/2 
X [ 7keV> 



B, 



CM B 



liiGf + (B/(xGy 



where ^S(A4) is a dimensionless shape function that rises 
steeply above M ~ 2.5 and plateaus to 1 above A4 ~ 10. 



In all work pres ented we use a fidu cial value of £ e = 0.005, 
as suggested inlHoeft et ah (2008), and the same as used 



Skillman et al.|(|20lTf 7 

There are several important things to notice about this 
model, which will help guide our interpretations of the 
results throughout this paper. First, the emission scales 
linearly with the downstream electron density, and with 
the downstream temperature to the 3/2 power. Addi- 
tionally, in regions where the magnetic field is less than 
the equivalent magnetic field strength from the CMB en- 
ergy density, Bcmb, the emission scales with B 1+S l 2 , 
where s ~ 3 for most relic situations. 

2.3. Analysis Tools 

In this work we relied h eavily on the data analysis and 
visualization toolkit, yt (Turk et al. 2011), to produce 
the derived data products presented. Here we describe 
the tools used specifically in our analysis, and leave fur- 
ther description to the yt documentation Derived 
Quantities^ such as WeightedAverageQuantity and To- 
talQuantity, are used to calculate weighted averages and 
totals of fluid quantities. For example, we use Weighte- 
dAverageQuantity to calculate the average temperature, 
weighted by cell mass. To analyze properties such as the 
radio and X-rayemission from our simulations, we use 
Derived Fields[_jto define the functional form of our new 
fields, which is then calculated on a grid-by-grid basis as 
needed. To calculate distribution functions, either as a 
function of position or fluid quantity, we utilize 1-D Pro- 
files and 2-D Phase Plot^j By specifying a binning field, 
we are then able to calculate either the total of another 
quantity or the average (along with the standard devia- 
tion). These are used to create radial profiles as well as 
characterize quantities such as the average magnetic field 
strength as a function of density and temperature. We 
also take advantage of adaptive slices and projections. 
Slicefj sample the data at the highest resolution data 
available, and return an adaptive 2D image that can then 
be re-sampled into fixed resolution images. Similarly, we 
use weighted and unweighted projections^] of quantities 
to provide average or total quantities integrated along 
the line of sight. Again, these adaptive 2D data objects 
can then be re-sampled to create images at various res- 
olutions. We also utilize and extend off-axis projections 
for use in integrating the polarization vectors of radio 
emission, to be described further in Section Finally, 
we use the spectral frequency integrator to calculate the 
X -ray emission based o n th e Cloudy code, as was don e 
Skillman et all (120111) and| Hallman fc Je l temal (|2011|), 



and was described in detail in |Smith et al.| ( |20"0g[ ). 

2.4. Polarization 

In addition to calculating the synchrotron emission, in 
this paper we investigate the polarization fraction and 
position angles of the emission. In order to compare our 
simulations to observations, we have developed several 

9 http://yt-projcct.org/doc 

10 http: / / yt-project.org/ doc / analyzing/objects. html#:derived- 
quantities 

11 http: //yt-project.org/doc/analyzing/creating_derived jields.html 

12 http:/ /yt-project.org/doc/visualizing/plots.html#:d-profiles 

13 http:/ /yt-project.org/doc/visualizing/plots.html#;slices 

14 http://yt-project.Org/doc/visualizing/plots.html#projections 
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new tools, including the ability to calculate the polariza- 
tion properties of the radio emission. In this section, we 
describe how we calculate Stokes I, Q, and U parameters 
from any viewing a ngle o f our s i mulation . For a review 



of the se topics, see Burn (1966); Longair (19941; Heiles 



(2002). While previous analyses ot polarized emission 



were capable of viewing along the coordinate directions, 
to our knowledge, this is the first presentation of 
off-axis polarized radio emission from AMR sim- 
ulations^ This capability presents several challenges. 
Whereas the total radio emission is calculated as a direct 
sum of the emission multiplied with the path length, the 
calculation of the polarized emission requires simultane- 
ous integration of each polarized component along the 
line of sight due to their mixing through Faraday rota- 
tion. 

We have built this capability on top of the analysis 
package yt. We began with the "off-axis projection'!^] 
operation, which is an off-axis ray-casting mechanism. It 
operates by creating a fixed-resolution image plane for 
which each pixel is then integrated through the simula- 
tion volume. To do this correctly, first the AMR hier- 
archy is homogenized into single-resolution bricks that 
uniquely tile the domain. This ensures that only the 
highest resolution data is used for a given point in space. 
These bricks are ordered and traversed by the image 
plane. The result of this is that we are able to inte- 
grate along the line of sight through the AMR hierarchy 
sampling only the highest-resolution cells for that given 
point in space. 

We have furthermore modified this framework such 
that the RGB channels of the image act as the total 
emission, I, and polarized emission along the x, I Xl and 



Iy , axes . 



I x and I y can be thought of as the emission- 
weighted electric field. The details of this calculation can 
be found in Appendix A. We first create derived fields 
that correspond to the magnetic field projected onto the 
unit vectors tlx, v y and Un, where iT x and v y are defined 
with respect to east and north vectors defined by the 
viewing direction, vti. We label these magnetic fields as 
B x , By, and Bu, respectively. We then define the po- 
larization angle \ of the electric field as the angle made 
between B x and B y rotated by tt/2. Finally, we define 
a Faraday rotation field Acj) — 2.62 x 10~ 17 x A 2 n e JB||dl, 
where all variables are in cgs units. Using a similar no- 



tation to Otmianowska-Mazur et ah] (20091), we then in- 

1/7 



tegrate along the line of sight the 
using the following discrete step: 



E) cos(Acj)) 



and 



I y values 



■ h+l - 




dl 


Ix,i+1 




dl f p (v x 


Jy,i+l- 




dl fp (Vy 









T (' 




Jy-.i- 





— sin(A(f>) 
cos(A(f>) 

where A<p is the Faraday rotation, v* x and v* y are the 
image plane coordinate vectors, f p is the fractional po- 
larization of the synchrotron radiation of a given power- 
law slope of electrons, dl is the ray segment length be- 
tween the incoming and outgoing face of each cell, and 
B is the magnetic field. Once integrated through the 



15 |Hoeft et al.| 12008) studied the view dependence on the total 
radio emission 

16 http:/ /yt-project.org/doc/cookbook/index.html#:cookbook- 
off axis-projection 



volume, we are then able to create intensity, polariza- 
tion fraction, and polarization direction maps. This ca- 
pability is available to download using the changeset 
with hash fc3acb 747162 here: https : //bitbucket . org/ 
samskillman/yt- stokes. 

3. RADIO RELIC PROPERTIES 

In this section we describe the general properties of the 
simulated galaxy cluster. We begin by comparing the 
morphological similarities between our simulated cluster 
and several observed clusters. We then move on to de- 
scribe the other gas properties in an effort to constrain 
the properties of the radio-emitting plasma. Finally, we 
will look at the time evolution of these quantities in order 
to understand their coupling during the merger process. 

3.1. Simulated Radio & X-ray 

We begin by comparing the radio and X-ray emis- 
sion from our simulation wit h the radio r e lics p resent 
in A3376 (see Figure la. 



Bagchi et alT] p006|)) and 



CIZA J2242. 8+5301 (see Figure 1 in Ivan Weeren et al 
(20101). In this work we calculate the X-ray emission 
using Cloudy to integrate the emission from 0.5 — 12 keV 
assuming a metallicity of Z /Z so i ar = 0.3. The resulting 
1.4 GHz radio and 0.5—12 keV X-ray emission is overlaid 
in Figure [T] The X-ray is shown in color with a dynamic 
range of 100. The radio flux is calculated by placing the 
simulated cluster at a distance of WOMpc/h. We then 
mask the radio emission such that 10~ 3 — 10 1 mJy is vis- 
ible. The total integrated flux for the left and right relics 
are 9.67x 10 24 and 3.11 x l0 2i W/Hz at lAGHz, which is 
simila r to many of the ob served single and double radio 
relics ( |Feretti et al.|2012 |. 



There are a few specific details that we highlight here 
due to their similarities to many observed radio relics. 
First, we note that this appears as a double radio relic. 
These are relatively rare compared to their single-sided 
counterparts. This snapshot is following a major merger 
roughly 300 Myr after core passage. The primary cluster 
is moving to the lower- left, with the secondary moving 
primarily to the right. After core passage, a merger shock 
develops, moving both to the lower left and upper right, 
and as will be seen in later figures, aligns with strong 
jumps in temperature and density. The two relic features 
are aligned with the direction of the merger as well as the 
shape of the X-ray emission. This alignment of the X-ray 



2012) 



Skillman et al. (|2011 


), as well as many 


van Weeren et al.|201 


lb 


Bonafede et al. 



The second key feature to this simulated double relic 
is the apparent aspect ratio of the radio emission. The 
length of the left relic (if we connect the two pieces that 
are separated by a short distance) is more than 2 Mpc/h, 
while the right-moving relic is 1.5 — 2.0 Mpc/h. While 
these relics are quite elongated, they are very thin. In 
many regions it is at most 100 — 200 kpc/h wide, even 
in projection. We note here, however, that this width is 
most likely underestimated since we are not tracking the 
aged populations of electrons that would exist for some 
amount of time behind the shock front. 

Comparing to A3376 (Ba gchi et al.|2006" ), we see many 
striking resemblances, including the complex morphology 
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Fig. 1. — Simulated X-ray and radio emission. The X-ray in the central regions shows a dynamic range of 100. 
calculated by placing the simulated cluster at WOMpc/h, and masked to show a dynamic range of 10 4 . 



The radio emission is 



of the Eastern portion of the relic as well as the ring- 
like structure to the outline of the radio emission. This, 
at the very least, suggests that the radio-emitting elec- 
trons are indeed related to the shock structures formed 
in merging galaxy clu sters. We see a similar structure in 
CIZA J2242. 8+5301 ( [van Weeren et al.)2010[ ), where the 
elongation of the X-ray emission points m the direction of 
the merger, aligning with the double radio relic. We also 
note the resemblance here with our simulation in terms 
of the very thin region of radio emission along the relic. 
This suggests that the cooling times of the relativistic 
electrons must be short. 

Figure [2] shows the fundamental quantities such as 
density, temperature, Mach number, and magnetic field 
strength, along with observable quantities such as the 
1.4 GHz radio and temperature fluctuations in the CMB 
due to the thermal Sunyaev-Zel'dovich effect. The ra- 
dio flux assumes a distance to the simulated cluster of 
100 Mpc/h. Each panel shows the same field-of-view 



and depth of 4.0 Mpc/h at z=0. In all panels except 
for the radio and Mach number maps, we have over- 
plotted the radio emission to help guide the reader's eye 
in determining the location of the emission relative to 
the underlying plasma. The top left panel shows the 
density-weighted temperature. Here the correlation be- 
tween the radio emission and the temperature structure 
is very strong, as the outward moving shocks are heating 
the gas to several x 10 8 K. Note that the sharp edges in 
the temperature structure, as well as the maximum in 
the temperature distribution, occurs 1.5 — 2Mpc/h away 
from the center of the cluster. This highlights the unre- 
laxed nature of this merging cluster. 

The bottom-left panel shows the density-weighted den- 
sity, and has a similar structure to the X-ray image shown 
in Figure [l] as is expected since the X-ray emission is a 
strong function of gas density. We note that unlike the 
temperature, the density is strongly peaked towards the 
center of the cluster, though the unrelaxed nature is ev- 
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ident by the elongation along the merger axis. The top- 
right panel shows the kinetic energy flux-weighted Mach 
number. In addition, we masked out all pixels that con- 
tributed < 10~ 10 of the peak radio emission. This was 
done to eliminate some of the shocks along the line of 
sight that are external to the cluster. The bottom-right 
panel shows a projection of the absolute magnitude of 
the magnetic field, weighted by density. Here we weight 
by density since in many regions there is no relic emis- 
sion, which would lead to difficult-to-interpret values in 
the those regions. Notice that the magnetic field also 
peaks in the center of the cluster around 1 /iG. There 
appears to be a drop in field strength outside the ra- 
dio relics. This makes sense because the merger shock 
is just reaching those regions, and then compresses the 
field behind the shock. 

We now use these four quantities to produce the images 
in the middle bottom column. First, the middle-bottom 
image s how s the radio emission calculated as outlined in 
Section |2. 2 Notice that as expected, the radio emission 



traces the regions that combine all four of the quantities 
in the left and right columns. This snapshot of the clus- 
ter properties suggests that radio emission requires a 
combination of dense, hot, magnetically threaded 
gas in the presence of moderately strong shocks. 
There are regions in the Mach panel of shocks in the 
A4 = 5 — 7 range that only contribute a small amount 
of radio emission. This region, which can be seen by 
rotating the viewpoint (not shown in this presentation), 
happens to lie outside the cluster a bit further where the 
magnetic field, temperature, and density have dropped 
considerably. We also notice that the strongest regions of 
radio emission correspond to Mach numbers in the 3 — 5 
ra nge with temperatures above 10 8 K , similar to findings 



(2011) 



in |Skillman et al. 

Finally, in the top-middle panel, we have overlaid the 
radio emission on top of the thermal Sunyaev Zel'dovich 
effect (tSZE). We have taken out the frequency depen- 
dence in the f(x) function, and multiplied the tSZE 
Compton y-parameter by the temperature of the CMB in 
order to get units of mK. Finally we multiply by 2.0, the 
maximum of f(x), to get the maximum decrement value. 
Notice the very strong correlation with the jump in tSZE 
and the presence of radio emission. Also note that this 
image is shown in linear-scale, and the dynamic range in 
this image is < 20. This has interesting implications for 
ongoing and futur e high-resolution tSZE me asurements 



with MUSTANG (Mroczkowski et al.. 
(IPlagge et al.||2012|), and CCAT (|Rac 



2011}. CARMA 



ord et al. 2009 ) 



Since the tSZF is sensitive to the integral ot the pressure 
along the line of sight, it may be preferable to X-ray stud- 
ies of galaxy cluster shocks because of it's linear depen- 
dence on density and temperature instead of a roughly 
quadratic dependence on density and sub-linear depen- 
dence on temperature. This also minimizes the effects 
of gas clumping and increases sensitivity to low density 
gas. This advantage is amplified in the outer regions of 
galaxy clusters where many of these radio relic shocks 
are located. 

3.2. Density & Magnetic Field Strength Relationship 

In the next two sections we will describe the physical 
properties of the radio emitting plasma in order to com- 
pare to what has been found observationally. In the fol- 



lowing analysis, we use the inner-most nested refined re- 
gion of the simulation (32 Mpc/h) 3 and ignore the lower- 
resolution regions in the remainder of the simulation. We 
choose to study this entire sub- volume instead of regions 
defined by the virial radius because the cluster is under- 
going a major merger. Within this region, we bin various 
quantities such as radio emission, density, temperature, 
and magnetic field. 

Before we discuss the radio emission, we first describe 
the structure of the magnetic field. In Figure [3] we show 
the average magnetic field strength as a function of den- 
sity, along with the 1 — a standard deviation where the 
average is weighted by density. T his is similar to plots 
found in Dubois & Teyssier ( 2008 ) , though in their study 
they also followed gas cooling. Co ntrary to what ha s 
been assumed in prior studies such as Hoeft et al. (2008); 
Skillman et al. (20111, we find that the magnetic held 
does not scale with p' 2/3 , but instead with p at low density 
and has an inner core with roughly flat magnet ic field at 
high d ensity. This was discussed previously in |Xu et ah] 
(20111, casting doubt on assumptions of any tight re 



lationship between density and magnetic field strength. 
This suggests that the manner in which we inject the 
magnetic fields may ultimately determine quantities such 
as the total radio luminosity. Future work should explore 
the effects of varying the magnetic field injection mecha- 
nism. For example, we may expect a shallower slope and 
more evenly distributed field strength if we inject mag- 
netic fields from multiple sources instead of seeding each 
cluster with a single source. 




Density [g/cm~3] 

Fig. 3. — Average magnetic field strength as a function of den- 
sity, where each cell is weighted by the density, with the shaded 
area denoting the standard deviation. For reference, we show an 
analytic function that is linear with p at low density and flattens 
at high density. 



Because of the steep drop in magnetic field strength 
at low densities, we expect that the radio emission will 
similarly decrease towards the outskirts of the cluster. 
Additionally, because the field strength flattens out at 
high density, we expect to see radio emission that is not 
necessarily biased to the highest density regions near the 
center of the cluster. 
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3.3. Kinetic Energy & Radio Emission Distributions 

The next quantity we will use to describe the gas prop- 
erties of the cluster are the kinetic energy flux through 
shocks and the radio emission based on the IHoeft fcl 
Briiggen (2007) model. We find it useful to show both 
of these quantities with respect to density, Mach num- 
ber, magnetic field strength, and temperature in order to 
characterize the gas that is most responsible for the con- 
version of kinetic energy to cosmic-ray electrons. For this 
particular study, we choose to show both kinetic energy 
flux as well as the radio emission. As seen in Section [272} 
the radio emission is a complex function of many quanti- 
ties, and it is difficult to disentangle each variable. The 
kinetic energy flux, however, is a much simpler quantity, 
relying only on the density, velocity, and temperature of 
the incoming gas: Fke = 0-5pv 2 Mc s . 

The top left panel of [4] shows the kinetic energy flux 
through shocks as a function of Mach number on the x- 
axis, and the magnetic field strength on the y-axis. Note 
that all un-shocked cells in the volume are ignored. For 
each bin, we calculate the kinetic energy, and normalize 
by the maximum value from all the bins. This distri- 
bution is shown in logarithmic space, as only a small 
number regions dominate the kinetic energy flux. In this 
figure there appears to be 3 primary populations of gas 
that contribute to the kinetic energy flux. The first is a 
very low Mach number (M < 1.5), high-magnetic field, 
region of the simulation, which corresponds to the turbu- 
lent flow in the center of the galaxy cluster. This slightly 
supersonic flow processes a large amount of kinetic en- 
ergy because of the high density and high temperature 
(and therefore sound speed) , even though the Mach num- 
bers are low. On the upper end of the shock Mach num- 
ber scale, there are shocks around a Mach number of 6-8 
with a magnetic field strength below 0.01 pG. These, 
as we will discuss later, are associated with shocks onto 
filaments. 

The third population in the kinetic energy flux is at 
Mach numbers between 3 — 5 with a magnetic field 
strength of between 0.1 — 1.0 /iG. These are the two 
primary merger shocks, which we will see are the main 
source of the radio emission. Note that this is not a sin- 
gle or even pair of distinct Mach numbers, meaning that 
observationally, a given merger shock should not 
be characterized by a single Mach number. If the 
properties of shock-accelerated electrons are strongly de- 
pendent on the Mach number, ascribing a single Mach 
number may lead to inconsistencies in the fitting to the 
observed radio emission. 

The lower-left panel shows the same quantity in color, 
but now decomposed into density and temperature bins. 
Here we get a different view of the same result. In this 
panel there are several knot-like regions in phase space, 
which likely correspond to each of the primary shocks 
in our cluster. However, in this case we see that there 
are 3-4 regions. There is a very low-density, 10 6 — 10 7 K 
region that likely corresponds to outer accretion shocks 
onto the filaments. The other regions of high kinetic en- 
ergy at more intermediate densities correspond to each 
of the primary merger shocks in the cluster. The large 
clump at very high temperatures and a density of roughly 
10~ 28 g/cm 3 corresponds to the left-moving shock at 
M. ~ 4. This will become even more apparent when 



we examine the right panels of this figure. The primary 
results of this study of the kinetic energy flux agrees well 
with prior stud ies of the distribution of kinetic energy 



flux in shocks dRyu et a l.||2003 Pfro m mer et al. 2006 
Skillman et al.||2UDHnVazza et al.|2009| |20TT] ). However 
this is the first time they have also been correlated with 
the magnetic field distribution. 

Now we contrast the distributions in kinetic energy flux 
with the radio emission. In the right panels of Figure |4j 
we use the color scale to signify the total radio emission. 
In the top right panel, we again show radio emission as 
a function of Mach number and magnetic field strength. 
Here we see that of the three features present in the ki- 
netic energy flux, only one remains in the radio emission. 
This corresponds to the regions where the three ingre- 
dients needed to create radio emission are all present. 
The Mach number is above the threshold for accelerat- 
ing high energy electrons, the gas is dense enough, and 
the magnetic field strength is high enough. We note that 
the falloff at low Mach numbers is due to the functional 
form of ^(M.) from Equation [l] For this simulation, we 
find that the majority of the radio emission in this 
cluster comes from shocks with Mach numbers of 
3.0—6.0 and magnetic field strengths of 0.1 — 1.0 /iG. 

The same general conclusions are found from the lower 
right panel. We see that, instead of a fairly broad dis- 
tribution in kinetic energy flux across densities and tem- 
peratures, the radio emission is confined to regions of 
relatively high densities and temperatures near 10 8 K. 
By combining this with the upper-right panel, we have 
determined the exact makeup of the gas responsible for 
the synchrotron radiating electrons in radio relics. 

It is important to note that the magnetic field in our 
simulation does not reach the levels calculated from ob- 
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van Weeren et al.||20l0| ). This suggests 
re magnetic field injection mechanism used in 
this study may not correspond to how it occurs in ob- 
served galaxy clusters. Alternatively, amplification of 
pre-existing fields prior to AGN injection, primordial 
magnetic fields, and several other processes such as tur- 
bulent dynamo (Iapichino & Briiggen 2012]) in t he post- 
shock region and the streaming instability (e.g. Achter- 
berg fc Wiersma|2007 1 may lead to higher magnetic field 



values, all of which could be investigated in future work. 
Finally, observations should be re-examined to determine 
if a lower value of magnetic field is possible, such as what 
has recently been found in Carretti et al. (2012). 



3.4. Merger Evolution 

Observationally, we are limited to a single snapshot in 
time, from which we must deduce the prior and future 
evolution of a given galaxy cluster. Fortunately this is 
not the case for simulations. In this section we describe in 
detail the evolution of the gas properties throughout the 
history of the clusters we are modeling. Figure [5] shows 
the evolution of bulk properties in the nested region of 
the simulation from redshift z = 2.95 to z = 0.0. There 
are five quantities shown here. 

1. Density: Average density, weighted by density. 

2. Temperature: Average temperature, weighted by 
density. 
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Fig. 4. — Phase plots of gas properties indicating the location of the kinetic energy flux and radio emission at shock fronts. The left plots 
show the kinetic energy distribution, while the right plots follow the radio emissivity. The top panels show the distributions as a function 
of magnetic field strength on the y-axis, and Mach number on the x-axis. The lower panels show them as a function of temperature on the 
y-axis, and density on the x-axis. 



3. Magnetic Field: Average magnetic field strength, 
weighted by density. 

4. Radio: Total 1.4 GHz Radio Emission 

5. Xray: Total 0.5-12 keV X-ray luminosity 

In each case we select the inner-nested region of the simu- 
lation, and use WeightedAverageQuantity or TotalQuan- 
tity "derived quantities" . Each of these averages and to- 
tals are then saved for future analysis. In the case of Fig- 
ure |5j we normalize each of the quantities by their maxi- 
mum value in order to fit them all on the same scale. In 
the case of the X-ray and radio luminosity, the emission 
is calculated using the blueshifted frequencies as these 
would be redshifted into the observer's frame to be at 
the correct values (0.5 — 12keV and 1.4GHz for X-ray 
and radio, respectively). 

We find a very interesting correlation with all of the 
fundamental and derived quantities. First, it is clear 
that there was not only the late-time merger near z = 0, 



but also earlier merger evolution near z = 1. First, we 
see that the density and temperature both start to rise 
0.2 — 0.5Gyr before the radio emission spikes. Analo- 
gously, the magnetic field and X-ray luminosity also fol- 
low this slow rise to a peak. Near the peak, the radio 
emission jumps up several orders of magnitude. This cor- 
responds to the formation of the shock front that then 
moves outwards from the cluster center towards the out- 
skirts of the cluster. After the core passage, the density 
and temperature returns to a lower but elevated level 
with respect to the pre-merger values. 

This suggests that the radio emission lags the merger 
event by a few x 10 8 years while the shock is setting up 
and expanding into the intracluster medium. It again 
highlights the dependence on not only the local charac- 
teristics of the emitting plasma, but also the shock sur- 
face area. Another key point is that the short timescales 
over which the radio luminosity varies implies that for a 
given mass or X-ray luminosity, there may be very large 
scatter in the radio luminosity. Overall, the radio emis- 
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Fig. 5. — Time evolution of integrated gas properties in the volume surrounding the structures of interest. For density, temperature, and 
magnetic field, we calculate a weighted average, using density as the weight. For radio and X-ray values, we calculate the total emission 
within the innermost nested region of the simulation. Note the simulation was run beyond z = to allow the merger to finish. 

sion is only above 10% of its peak value for ^0.5 Gyr. 
This could help explain the observed lack of radio emis- 
sion f rom clusters with o bvious mergers such as Abell 
2146 (Russell et al. 2011). One caveat to this result is 



that we do not follow the electrons as they cool. How- 
ever, because the cooling timescale, 

r « 2 x 10 12 7" 1 ((1 + zf + (B/S.SfiG) 2 )- 1 years (3) 



tore 



for these el ectrons with 7 ~ 1000 — 5000 is short |Ferrari 
et al. |2008|), this additional time has little effect. 



There 

e characteristics of this time evolution should not 



change substantially with a proper treatment of the ag- 
ing electron population. 

4. OBSERVATIONAL IMPLICATIONS 

In this section we set out to provide a theoretical per- 
spective on the analysis of observed radio relics. In par- 
ticular, we comment on some of the assumptions that are 
often employed, and in several cases point out how these 
may be dubious. We begin by examining how the view- 
ing angle of a radio relic can impact its interpretation. 
We expand on this point in the context of spectral index 
analysis, where not only does viewpoint play a role, but 
small-scale fluctuations in the shock properties can lead 



to observational signatures that mimic an aging popula- 
tion of cosmic ray electrons. Finally, we point out the 
limitations of our models; specifically, we must include 
aging populations of electrons in future calculations in 
order to capture accurate polarization fraction and di- 
rection. 

4.1. Viewpoint 

Observationally, we are limited to a single viewing an- 
gle for each object. Unfortunately, because radio relic 
emission is not spherically symmetric, there will be a 
viewing angle-dependent emission strength and morphol- 
ogy. In this section we set out to demonstrate this fact 
and how it can affect our interpretation of radio emit- 
ting regions in galaxy clusters. We begin by taking our 
original viewpoint and rotate the viewing angle by 180 
degrees over 18 frames. The result is shown in Figure 
[6j where radio emission is projected along each viewing 
angle. 

What can be seen from Figure [6] is that, while for some 
orientations the radio emission forms an obvious "double 
relic" configuration, other orientations yield what seems 
to be a single, more diffuse, object. However, given a 
single frame, it would be difficult to determine the true 
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structure of the emission. In fact, this may lead to a mis- 
interpretation of the radio emission to be that of radio 
halo origin . W e will expand upon this line of reasoning 
in Section |4.2| There are several other things to note 
here in the rotation of the radio emission. We see that 
in some orientations (see top row, 5th & 6th columns), 
we reproduce morphologies that include two outer relics, 
with one of the relics ending up between the two, close to 
the center of the cluster. This is very si milar to observed 
clusters such as MACS J1752.0 + 4440 dvan Weeren et al.| 



2012 1, CIZAJ2242.8 + 5301 (van Weeren et al. 2 011d[ ), 
ancTMACS J0717.5 + 3745 ( |van Weeren et al.||200y| . It 
may be possible that the emission in these clusters be 
not of radio halo origin, but simply radio relic emission 
viewed coincident with the cluster center. Similar work 
has been done with hydrodynamic simulations of galaxy 
clusters and viewi ng them along the coordinate axes in 
Vazza et al. (2012), where they found that the low emis- 
sivity due to the small size along the line of sight may 
explain the lack of central radio relics. 

4.2. Spectral Index 

Observationally, the spectral index of cluster radio 
relics is measured by comparing the emission at several 
different observed frequencies. In our work, we calculate 
the spectral index by first calculating the radio emission 
at two frequencies (here we use 1.4 GHz and 330 MHz). 
After smoothing by a Gaussian with a full-width half 
maximum size of 4 pixels, we then use these two maps 
to calculate the spectral index, similar to what would 
be done observ ationally. Because we use the |Hoeft fc| 



Briiggen (20071 model, we will recover a spectral index 
from comparing these two maps that is equal to the cu- 
mulative spectral index due to the emission from elec- 
trons over their entire lifetime. This is steeper than 
the spectral index that would be predicted from the 
prompt emission from a shock front, which is related by 

^prompt — (1 ^Qintegrated) /2- 

In reality, what is thought to happen is the leading edge 
of the shock front should accelerate electrons to aiprompt ■ 
As the radio emitting electrons move downstream from 
the shock they cool and the spectrum steepens. Qual- 
itatively, we would expect that "edge-on" observations 
should produce values near a promp t, whereas a "face-on" 
view would be sensitive to the entire lifetime of the elec- 
trons, and therefore be closer to on n t e grated- 

Spectral stee pening is found to happe n in several ob- 
served clusters ( van Weeren et aL]|2011a ) and measuring 
the steepening of the electrons as they progress away 
from the shock can help constrain the loca l magnetic 
field, as was done in |van Weeren et"aL (2011a I. However, 
the assumption in these calculations is that the spectral 
steepening is due entirely to the aging of the electrons. 
However, we see similar steepening even though we do 
not include the spectral aging of the electrons as they 
advect downstream! These variations are entirely due 
to a varying magnetic field and shock strength in a non- 
uniform medium. The spectral indices of the prompt and 
integrated spectra are shown in Figure[7j and the fluctua- 
tions in the underlying fields for the lower-left "edge-on" 
relic are shown in the right panel of Figure [8j 

In Figure [7j we see that the "edge-on" view shows a 
spectral steepening whereas the "face-on" view shows the 
spatially variant spectral shape due to the shock proper- 



ties. For completeness we have mapped the colormap 
of both viewing angles to both the prompt and inte- 
grated spectral indices. The same qualitative behavior is 
seen with both calculations. These fluctuations are ex- 
plained in Figure [8j where we see that the radio emission- 
weighted Mach number can vary between 2 — 8, and the 
magnetic field can vary from 0.1 — l/.iG. Therefore we 
warn that extrapolating from the measured spec- 
tral steepening to calculate gas properties should 
be done with care, as projection effects and spa- 
tial variation of the gas properties are important. 
It may be insufficient to prescribe a single Mach 
number or magnetic field to a given radio relic. 

4.3. Polarization Fraction & Direction 

Using our newly developed method, we have created 
radio polarization fraction and direction maps for our 
galaxy cluster from two viewing positions. The first gives 
an "edge-on" view of the radio relic, whereas the second 
is aligned such that the relics are viewed "face-on." We 
will examine, in detail, the polarization signatures of the 
"edge-on," and then describe the differences in the "face- 
on" view. 

In Figure [9j we show the linear polarization fraction 
and direction of 1.4GHz radio emission. Lines denote 
the local linear polarization direction, the length of which 
corresponds to the polarization fraction. For clarity, the 
polarization fraction is also shown in color. The polar- 
ization fraction is calculated using 



fp = 



1 



(4) 



In this case we use a 512 x 512 pixel image plane to 
project through a cube of length 4 Mpc/h on a side. 
At full resolution (7.8kpc/h), we see that the polariza- 
tion fraction reaches a maximum of ~ 75%, and that the 
polarization direction is correlated along the relic, pri- 
marily perpendicular to the shock. This is very similar 
to what is found observat ionally in CIZA J2242. 8+5301 
(|van Weeren et aL"||2011a ). The authors find strong (50- 
757o) polarization at what is presumed to be the leading 
edge of the shock. Additionally, they find that the polar- 
ization direction is fairly constant over the length of the 
relic. We do, however, see greater variation in the po- 
larization fraction and direction both across and along 
the relic. These fluctuations in our simulation may sug- 
gest that there are additional physical processes which 
lead to a more ordered field. In order to investigate the 
small scale fluctuations in the polarization direction, we 
examined a slice of the magnetic field strength and direc- 
tion through this relic. This is shown in Figure [8} What 
we found is that while the shock (overlaid in white) cuts 
through regions which have fairly strong variations in the 
magnetic field direction, the region "behind" the shock 
has a magnetic field that is compressed along the shock 
propagation direction. Therefore, it may be possible that 
if we were to follow the evolution of the cooling electrons 
as they moved downstream across the shock, the mag- 
netic field that they reside in may become more ordered 
parallel to the shock, leading to a longer coherence length 
as more constant polarization direction. Therefore we 
would expect that in future work when we exam- 
ine the emission from an evolving population of 
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Fig. 7. — Spectral index of simulated radio relic emission. The left portion of the image shows the "edge-on" view, whereas the right 
shows the "face-on" view. Both views are on the same scale. The left colorbar shows the mapping of color to the integrated spectral index 
including particle aging. The right colorbar shows the mapping of color to the prompt spectral index. Both colorbars apply to both views, 
providing a rough estimate of the uncertainty in our models of the spectral index. 
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Fig. 8. — A zoom in of the lower left radio relic. The left two panels show radio emission-weighted projections of the Mach number (left), 
and magnetic field strength (middle). The right panel shows a slice of the magnetic field strength, with black lines indicating the local 
magnetic field direction in the plane of the slice and the white overlay show the location of cells identified to be shocks. 



cosmic rays, we should see even better agreement 
with observations of polarization direction. 

In the top-right panel of Figure j9j we show the same 
polarization map, but this time for when the relic is 
viewed "face-on" . We again see a very high polarization 
fraction. However, this time the polarization direction is 
significantly less coherent. This is due to the magnetic 
field not being modified as strongly in the plane of the 
shock as it is perpendicular to the shock. Therefore the 
turbulent structure is preserved in the image plane and 
the polarization vectors are not preferentially modified. 

However, the behavior of the polarization fraction and 
direction drastically changes if we then apply a guassian 
kernel with a size of 4 pixels. In the "edge-on" view, 
the polarization fraction and direction is fairly well pre- 
served. The fraction only drops to between 30 — 65% 
and the direction is still fairly correlated across the relic. 



In contrast, for the "face-on" view the polarization has 
dropped to between — 15% in most regions. This is a 
classic example of beam depolarization. Because the po- 
larization direction is highly disordered, smoothing the 
image drastically reduces the overall polarization. 

Currently our simulations do not exhibit as much of 
a constant polarization as that found in observations of 
clusters such as CIZA J2242. 8+5301, where the polar- 
ization direction is constant over Mpc-scale distances. 
There are several possible explanations. Because we are 
not tracking the electron distribution as it cools behind 
the shock, we may be missing the emission from the more 
ordered field line regions behind the shock. Simulations 
capable of tracking these electrons are therefore needed 
to explore that possibility, and will be addressed in future 
work. If doing so is still incapable of producing Mpc-scale 
ordered polarization maps, it may suggest that the injec- 
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Fig. 9. — Polarization fraction and direction. In each panel, the polarization direction is denoted by the black quivers, while the 
polarization fraction is represented by both the color scale as well as quiver length. The top panel shows the relic at full resolution 
(7.8kpc/h), while the lower panels show the same view at 4 times worse resolution. At z = 0.2, the redshift of CIZA J2242. 8+5301, this 
corresponds to angular resolutions of 3.36" and 13.44", respectively. The left panels show the polarization for the "edge-on" (top) while 
the right shows the "face-on" view. 



tion mechanism or magnetic field evolution is different 
than what we have simulated. 

5. DISCUSSION & FUTURE DIRECTIONS 

We have carried out high resolution MHD AMR cosmo- 
logical simulations using an accurate shock finding algo- 
rithm with a radio emission model for shock-accelerated 
electrons to examine the properties of radio relics in 
galaxy clusters. We summarize the physical conditions 
of the cluster, and several of the warnings for the inter- 
pretations of observed radio features: 

• Cosmological initial conditions lead naturally to 
the formation of giant radio relics whose proper- 
ties are very similar to observed relics. This is a 



natural result of the process of mergers that create 
galaxy clusters, where the magnetized intracluster 
medium is subject to a series of strong, large-area 
shocks. 

• We find that the radio emission in our simulated 
clusters is strongest in plasma that has 0.1 — 1.0 /iG 
magnetic fields with a shock Mach number between 
3 — 6, with densities near 10 -28 g/cm 3 and temper- 
atures of 10 8 K. 

• If the shock acceleration efficiency for electrons is 
higher than assumed at low Mach numbers, a large 
reservoir of kinetic energy becomes available in hot 
(10 7 K) gas at both lower and higher densities than 
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what is currently producing emission in our simu- 
lations. 

• We find a magnetic field distribution that does 
not follow the often assumed B oc p 2 / 3 relation. 
Instead, it follows B oc p at low densities, and 
flattens out to sub — fj,G levels above densities of 
l(T 27 g/cm 3 . 

• Warning: Without a model of spectral aging of 
electrons, we still recover a spectral index gradi- 
ent, indicating that observed gradients should not 
necessarily be interpreted as the spectral aging of 
electrons. It may simply be due to the projection 
of a curved shock front with varying Mach number. 

• Warning: From the time evolution of our simu- 
lated galaxy cluster, it is clear that only a small 
portion of its lifetime may be spent in a regime 
where it is bright in the radio wavelengths. This 
may explain the apparent lack of observed radio 
emission in some massive galaxy clusters with early 
or late phase merger. 

• Warning: The viewing angle of a given galaxy 
cluster may have significant impact on the classi- 
fication of its radio emission. Double radio relics 
viewed at some orientations are difficult to differ- 
entiate from radio halo emission. 

• Warning: Future simulations of galaxy cluster ra- 
dio relics must follow the temporal evolution of 
the electron population in order to reproduce valid 
polarization results, as the downstream conditions 
from the shock may include more ordered magnetic 
fields, leading to different polarization maps than 
when the emission is assumed to come only from 
the shock center. 

Many of these warnings apply both to observational 
and theoretical studies, and the classification and anal- 
ysis of radio relics in all contexts must be done with 
care to avoid confusion with radio halo emission. There 
are several advancements that can be made theoretically. 
We are in the process of developing and testing the nu- 
merical framework necessary to follow the cosmic ray 
electr on and proton popu lations, using a meth od simi- 
lar to piniatj] ([200l]) and [JoneTfe^Ka^ ([2005]) . Doing 
so will allow us to probe the spectral distribution of these 
non-thermal populations in the context of high- resolution 
cosmological simulations. Once this is merged with our 
ability to produce synthetic emission and polarization 
maps, we will be able to more directly compare to cur- 
rent observations. Additionally, we do not explicitly in- 
clude any magnetic field source terms at the shock fronts. 



Exploring local fiel d generation mechanisms such as the 



Weibel instability (Weibel 1959 Fried 19591 may allow 
for alternate magnetic field strength and structure be- 
hind the shock front from what we present here. Finally, 
we may need to push to higher resolution studies at the 
shock fronts to be able to follow the shock-amplification 
of magnetic fields. Doing so in a cosmological simulation 
is currently intractable; however, improvements in com- 
putational speed and parallelization may allow for future 
studies. 

Observationally we are entering a golden age of ra- 
dio telescopes with the upgraded Jansky VLA, GMRT, 
and LOFAR, and are looking forward t o the SKA, and 
possibly a lunar farside radio telescope! Burns & Lazio 



2012 Lazio et al. 2009). These improvements will lead 



to greater sensitivity and bandwidth, allowing for multi- 
frequency studies of galaxy cluster environments in un- 
precedented detail. However, only through a coordinated 
effort between simulations and observations will we be 
able to fully understand the plasma physics in these cos- 
mic environments. 



thank the referee for in- 
a much improved paper. 
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APPENDIX 
POLARIZED EMISSION INTEGRATION 



In this section we detail the integration methods used to c alculate the tot al an d polarized radio emission Faraday 
rotation. Where possible, we follow the conventions listed in Longair (1994) and |Otmianowska-Mazur et al. (2009). 

As input to the system, we will assume that there is an emissivity, e [erg /s/ cm 3 /Hz}; magnetic field, B [67]; electron 

number density n e [cm~ 3 ]; and viewing angle, L. 
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From the viewing angle and magnetic field, we first decompose the magnetic field into components along the line of 
sight Si I and perpendicular to it B±. We use 6 to represent the angle between the L and B, where we 9 has a range 
of0-7r. 



S|| = B-L = \B\cos(9) 
B ± = \B\sin(9) 



(Al) 
(A2) 



We then further decompose B± into components aligned with the east and north vectors in the image plane, here 
referred to as B x and B y . Using B x and B y , we then define an angle field, Xi 



X = (tan-^By/Bx) + tt/ A) % tt 



(A3) 



where we explicitly cast the tan 1 into [0..7r]. \ represents the local direction of the electric field in the image plane. 
We then use this value to calculate the fractional polarized emission along the east and north vectors, 

(A4) 
(A5) 

(A6) 



/ p = (« + l)/(« + 7/3) 
t P x = ef p \ — \cos(x) 



where f p is the polarization fraction, defined using the spectral index of the relativistic electrons, s. 

Once these emission terms have been defined, we then integrate through a simulation, accounting for the Faraday 
rotation, 



Aip = 2.62 x 10" 



'n e \ z B\\dl 



(A7) 



where dl is in units of cm. All of this is then integrated to calculate the total intensity /, and the polarized components 
I x and I y , which is written using the discrete approximation: 



h+l 




Ix,i+l 









dl 

dl f P (vx • E) cos(A<p) —sin(A(j)) 
dl f p (vy ■ E) sin(A<p) cos(A<j)) 



(A8) 



where the index i represents the i'th cell along the line of sight. 

As a test of this method, as implemented using yt, we initialize a completely polarized background at the simulation 
domain boundary, with the polarization angle pointing horizontal in the image plane. To do so, we override the 
emission, polarization fraction, density, and magnetic field strengths in the simulation presented in this paper. In this 
case, we set the magnetic field equal to (though for numerical reasons we use a very small number), except for a 
sphere of with a radius equal to one-quarter of the simulation domain size. Inside this sphere, we allow the magnetic 
field parallel to the line of sight to be exactly that which will, for rays passing through the center of the sphere towards 
the "observer", be enough to rotate the polarization vector by tt. We show the results of this test using 128 2 pixels 
in Figure 10 for both on-axis (L = (1,0,0)) and off-axis (L = (1.0,0.3,0.7)), where we find excellent agreement with 
the expected behavior. For both viewing angles, the maximum of f p — 1.0 ~ 10 -14 , and the polarization angle has 
a fractional error in the center of the image of 1.7 x 10~ 5 for on-axis, and 1.3 x 10~ 2 for off-axis. Given the coarse 
nature of the sampling of this image and the likelihood that the center pixel doesn't exactly traverse the center of the 
sphere, these should be viewed as upper limits. 

As a second test, we initialize two y — z planes of radiation each with electric field vectors that are perpendicular 
to each other along. For the first plane, we set B x = B z = 0.0, and B y = 1.0. The second plane has B x = By = 0.0 
and B z = 1.0. Each plane has a width equal to ^ of the simulation width, centered at x = and || in units of 
the simulation width. We set the image width equal to 2.0 and use 256 2 pixels. Finally, we begin by looking from 
(—1,0,0) towards (1,0,0), centered on (0.5,0.5,0.5), and rotate in the x — y plane. Therefore in all cases plane 2 
has a magnetic field that is vertical and plane 1 begins with a horizontal magnetic field but then proceeds to have 
a component along the line of sight. For this test we disable the Faraday rotation, and examine only the effects of 
combining two planes of radiation on off- axi s lines of sight. We vary the rotation from the original viewpoint to have 
(0,5,15,30,60). 



angles of 9 



In Figure 



11 



we show the fractional polarization as a function of image pixel. In all 
cases for the portions of the planes that overlap along the line of sight, the calculated polarization fraction is within 
double precision fractional error of the expected polarization fraction. 

All of these capabilities are demonstrated in a public repository, found using the changeset with hash fc3acb747162 
here: https://bitbucket.org/samskillman/yt-stokes. Future improvements to this code as well as tighter inte- 
gration in the primary yt repository is expected, both in performance and usability. 
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Fig. 10. — On and off-axis Faraday rotation test. The left panel shows the on-axis Faraday rotation through a magnetized sphere, with 
an image width equal to the domain size. The right panel shows the same rotation, but off-axis and with a width of 1.6 larger than the 
left, to show the off-axis nature of the domain. The electron number density and magnetic field strength of the sphere are chosen to rotate 
the polarization angle 7r radians for the rays passing through the center of the sphere. 




1.5 

Image x pixel 



Fig. 11. — Dual plane polarization test. The polarization fraction as a function of image pixel across the mid-plane of the image, shown 
for varying viewing angles that are measured as an offset in the x — y plane from a viewing direction of L = (1, 0, 0). 



